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ABSTRACT 


N  The  problem  considered  is  that  of  finding  the  thermal 
entrance  or  thermal  mixing  region  of  a  pipe  where  the  wall  heat 
flux  is  zero  along  the  length  of  the  pipe.  The  fluid  is  assumed 
to  enter  the  pipe  with  a  non-uniform  temperature  profile  and  a 
fully  developed  turbulent  velocity  profile. 

The  approach  is  analogous  to  that  introduced  by  Pettis  for 
the  solution  of  the  uniform  wall  temperature  problem. 

The  results  for  different  radii  tubes  and  Reynolds  numbers 
are  presented  in  graphs  which  show  the  adjustment  of  the  temper¬ 
ature  profile  down  the  pipe.  These  results  provide  information 
concerning  the  minimum  length  mixing  tube,  since  any  wall  losses 
will  require  a  greater  length. 
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n 


exp 
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k 


Pr 
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arbitrary  constants 
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specific  heat 
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pipe  diameter 
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parameter  defined  by  Equation  27 


total  heat  flux 
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flux  along  the  inner  wall  surface 
pipe  radius 

Non-d Imens lonal iz ing  pipe  radius 

radial  coordinate 

non-dimensional  radius  r 
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constants  defined  directly  under  Equation  25 
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decay  factor  for  the  temperature  in  the  arlal  directijn 
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indicates  partial  differentiation 

temperature  difference  betveen  the  local  temperature 
and  the  temperature  at  the  inner  vail 

temperature  at  x  *  0  and  r  *  0 

entrance  temperature  distribution 

eigenvalues  for  Equation  13 

dynamic  viscosity 

kinematic  viscosity 

density  of  the  gas 

refers  to  a  summation 
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I.  INTRODUCTION 

The  problem  of  heat  transfer  in  fully  established  turbu¬ 
lent  flov  in  cylindrical  tubes  has  received  considerable 
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attention.  In  all  cases  the  flov  Reynolds  number  is 

sufficiently  large  to  Justify  the  assumption  of  negligible  axial 
conduction  in  the  fluid.  As  a  consequence,  many  of  the  mathe¬ 
matical  investigations  were  reduced  to  solving  an  eigenvalue 
problem  vith  the  vail  temperature  taking  the  form  of  a  step 
function.  Once  this  vas  accomplished,  other  boundary  conditions 
such  as  prescribed  heat  flux  or  prescribed  vail  temperature 
variation  vere  included  by  the  method  of  superposition.  The 
case  of  a  perfectly  insulated  tube  vail  is  considered  in  this 
paper . 

The  fluid  enters  the  pipe  vith  a  non-uniform  temperature 
and  a  fully  developed  turbulent  velocity  profile.  At  succeed¬ 
ing  axial  stations,  due  to  the  radial  conduction  and  turbulent 
mixing,  the  temperature  vill  deviate  from  the  entering  profile 
until  a  uniform  distribution  is  approximately  achieved.  Pipes 
vith  this  length  are  used  as  mixing  tubes  and  the  adjustment 
or  damping  of  the  temperature  profile  is  important  in  estimating 
the  exit  temperature  profiles  of  shorter  tubes. 

The  analysis  is  similar  in  the  general  mathematical  approach 
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to  that  presented  by  Latzko  and  Pettis  for  an  isothermal  vail. 
They  assumed  a  one-seventh  pover  velocity  profile,  a  simplified 
eddy  diffusivity  and  a  Prandtl  number  of  unity.  Hovever, 
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Latzko1  obtained  crude  approximations  to  the  first  three 

p 

eigenvalues  by  using  Legendre  polynomials,  vhile  Fettis  obtai 
ed  good  estimates  of  the  same  three  eigenvalues  by  the  use  of 
Jacobi  polynomials. 
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II.  STATEMENT  OF  THE  PROBLEM 


A  schematic  diagram  shoving  the  coordinate  system  is  given 
in  Figure  1.  We  shall  consider  the  section  of  pipe  to  the 
right  of  x  *  0,  where  the  vail  heat  flux  is  equal  to  zero.  The 
flow  possesses  a  fully  developed  turbulent  velocity  profile 
and  a  selected  entrance  temperature  profile  at  x  ■  0. 

Subject  to  the  limitations  given  below,  the  steady  state 
energy  equation  is: 


.  dO  I  d.  % 

r  aTrH*  (i) 

vhere  qr  is  the  radial  heat  flux,  positive  in  the  +r  direction. 

In  writing  the  energy  equation  (l),  the  usual  basic  assump¬ 
tions  are  adopted: 

a.  The  fluid  properties  are  assumed  constant. 

b.  Viscous  dissipation  is  negligible. 

c.  Axial  diffusion  of  heat  is  negligible  compared  to  the 
axial  convection. 

d.  The  flow  is  hydrodynamlcally  fully  developed. 

The  statement  of  the  problem  is  completed  when  the  initial 
and  boundary  conditions  are  specified  for  the  function  0(r,x). 
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At  the  inner  surface  of  the  pipe,  ve  have: 


where  kv  and  q^  represent  the  gas  conductivity  and  flux  at 
the  wall,  respectively. 

The  fluid  temperature  at  the  entrance  (i.e.,  at  x  *  O) 
is  given  by  the  following: 


(2d) 


Equations  1,  2  and  2a,  and  also  the  imposed  condition 
that  no  infinite  temperature  exists,  constitute  the  mathematical 
statement  of  the  problem. 
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III.  AHALYSIS 


In  a  previous  paper,  Latzko's  differential  equation 
from  Equation  (1) 


for  convective  heat  transfer  in  fully  developed  turbulent  flov 
vith  a  velocity  profile  represented  by  the  equation 


(3a) 


vas  solved  subject  to  the  following  boundary  condition  at  the 
vail: 


qr*Ur0 


The  present  paper  considers  the  same  equation  (3),  but  vith  the 
boundary  condition  q  *  0  at  the  vail.  This  condition  when  the 
inlet  temperature  profile  is  constant  leads  to  a  trivial 
solution.  In  order  to  obtain  a  non-trivial  solution,  an  initial 
temperature  represented  by  the  function 


©«  a,  (i  -e,*14 
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1 8  used,  vhere 
and 


F  Is  defined  at;  the  non-dimensional  radius.  — , 

R 


Z  ■ 


A  solution  to  Equation  3  can  be  vritten  in  the  form 


0*g(r)  exp(-£x) 

and,  since  the  fluid  temperature  and  vail  temperature 
asymptotically  approach  each  other  as  the  pipe  length  increases, 

0  must  be  positive.  When  the  solution  as  given  by  Equation  6, 
is  inserted  into  Equation  3  and  use  is  made  of  the  non-dimension¬ 
al  radial  transformation  previously  defined  as  z,  the  differen¬ 
tial  equation  (3)  becomes 


jL 

dz 


*-wz7g 


(7) 


vhere  oj  is  defined  as 


p 

A  procedure  analogous  to  that  used  by  Fettls  for  the 
Isothermal  vail  problem  vas  folloved  in  obtaining  a  solution  to 
equation  7  for  this  case,  vhere  the  surface  heat  flux  is  zero. 
The  appropriate  boundary  condition  is: 


ag 

dz 


*0  at  z*0 


(8) 
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In  addition,  because  of  the  singularity  in  Equation  7,  i.e., 
vhen  z*l,  it  is  necessary  to  require  that  the  function  g  be 
finite  at  the  centerline  of  the  pipe. 

Solutions  of  Equation  7  which  satisfy  the  above  requirement 
and  boundary  condition  will  exist  only  for  discrete  values  of  a 
and  have  the  form 


0* 


^9nWexp(-^)x 


(9) 


where  g(z)  is  the  eigenfunction  solution  of  Equation  7.  The 
constants  "D"  must  be  determined  to  satisfy  the  initial  condition 
0(2,o)*0Oi  where  Qq  is  prescribed  entrance  temperature  profile. 
This  may  be  accomplished  by  use  of  the  orthogonality  property  of 
the  eigenfunctions,  gn(z),  that  is 


«nlz,dz  *° 


when  m  t  n«  Setting, 


ZDn«n<*>*®o 


n*o 


(10) 


(ID 
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ve  find  that 


D  - 

/o'l7«nt2)  dz 


(12) 


Before  obtaining  the  solution  of  Equation  7 ,  ve  vish  to 
consider  the  solution  of  the  auxiliary  equation 

<l3> 

vhich  is  similar,  but  not  identical,  to  Equation  7.  With  the 
change  of  variable  t*z 7  ,  Equation  13  becomes  an  equation  of 
the  hypergeometrlc  type: 

(,4) 

The  details  for  obtaining  the  eigenvalues  and  eigenfunctions 
for  Equation  lU  are  given  in  the  Appendix. 

Nov  assume  that  a  solution  of  Equation  7  can  be  found  in  the 

form 


gn<z>'£)>n  fcnte) 
n=o 
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If  Equation  15  is  substituted  into  Equation  7,  ve  obtain 


Because  the  eigenfunctions  $n  also  satisfy  the  auxiliary 
Equation  13 »  Equation  16  becomes 

^OnfxnzB®n(z)-wnz7<X'n{j)l*0  (17) 

n=o  L  J 

We  can  now  make  use  of  the  Galerkin  method  to  determine  the 
Dn  by  requiring  that  the  left  side  of  Equation  17  be  orthogonal 
to  the  $n  for  n  *  0,  1,  2  ...N,  thus  arriving  at  the  following 
system  of  equations: 


for  m  *  0,  1...N,  with  Aan  defined  by  the  equation 


(18) 


Amn  «jfz7$m(z)<!>nU)dz 


(19) 
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The  characteristic  equation  of  the  system  given  by  Equation  18 
is  the  determinant 


I  ^mn ~ ^mnl  *  0  (20) 

n=o  *■  -* 

The  roots  of  Equation  20  give  approximations  to  the  first  n 
eigenvalues  of  Equation  7,  and  the  complete  solution  to  Equation 
3  is  given  by  Equation  9* 

In  the  present  case,  where  the  •n<  z)  and  therefore  the 
gn(z)  are  expressed  as  polynomials,  the  coefficients  Dn  can 
be  obtained  explicitly  provided  the  initial  temperature  distri¬ 
bution  can  also  be  described  by  a  polynomial.  For  example,  a 
suitable  entrance  temperature  profile  could  be  represented  by  the 
function 


6jf(r)=0j  (l-r2)2=ffjZ14 


(21) 


By  equating  like  powers  of  the  function  z7  ,  we  obtain  the  follow¬ 
ing  system  of  algebraic  equations  for  the  Dn : 


+0.26688  D, 
-0.31298  D, 
-0.86258  D, 
+  0.51278  D, 
-0.18864  D, 


+  0.22691  D2  +  0.19413  D3  +  0.32980  D4  *  0 

-  0.68164  D2  -  0.98781  Dg  -  5.42133  D4  *  0 

-  1.51668  D2  -  2.83335  Dg  +  21.29747  D4  »  1 

+  2.63729  D2  +  11.27409  Dg  -  30.31206  D4  =  0 

-  0.36071  D2  -  7.95572  Dg  +  14.22400  D4  =  0 


(22) 
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Solving  the 

last  four 

of  Equations  22  for  the  coefficients  Dn , 

ve  obtain 

I 

D,  = 

-1.19346 

1 

* 

d2  = 

0.25817 

•  (23) 

* 

0.061657 

d4  * 

0.025205 

Substituting  these  four  values  for  the  coefficients  Dj  through 
D4  into  the  first  of  Equations  22,  we  obtain  for  DQ  the 
numerical  value 

o 

Dc  «  0.23965  (24) 

The  complete  solution  as  given  by  Equation  9  can  now  be  vritten 
in  the  form 

g-'Do+  **p(fif)+D2  92 

('  |r  *)  +  DJ  g5  e‘p(*  pftx) 

+  04«4  eXp(|lX) 


(25) 
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^qT 


or 


*0.23965-1.9346  exp  (- -|Lx)|o.26688-0.3I298z7 

-0.86258Z14  +  0.5I278Z21  -O.I8864z28] 

,2 


22691  -  0.68I64Z7 


+  0.25817  exp 

'  #  l 

-  1. 516682 14  +  2.63729z*‘  -0.360712**1 

+  0.06166  exp  (•>)[»  .19413 -0.9878IZ7 

-2.83335Z14  +  II.27409Z*1  -  7.95572z“l 

+  0.025205  exp  (-  -|g-  x)[o.32980  -  5.42I33Z7 

+  2I.29747Z14  -30.3I206Z21  +  I4.22400Z**  j 


(26) 


where 


P* 


8x2  $ 
7x0.199 


(27) 


and  F>  represents  the  inner  wall  radius. 
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IV.  DISCUSSION  AND  CONCLUSIONS 


The  approach  employed  in  obtaining  a  solution  to  Equation  3 

for  the  specified  boundary  condition  was  analogous  to  that  first 
2 

used  by  Fettis  in  solving  the  isothermal  problem  with  a  uniform 
entrance  temperature  profile.  However,  to  obtain  a  non-trivial 
solution  for  the  present  case  it  was  necessary  to  assume  that 
the  fluid  entering  the  pipe  had  a  non-uniform  temperature 
profile.  The  eigenvalues  a»n  and  the  eigenfunctions  ^  remain 
unchanged  for  any  profile  chosen,  because  neither  depend  on  the 
initial  temperature  distribution.  In  the  chosen  numerical  example, 
the  initial  temperature  distribution  was  represented  by  the 
following  equation: 


0(F,o)=0,f(r)=0j(|-r2)2  (28) 

where  0j  is  the  initial  centerline  temperature.  The  centerline 
temperature  as  a  function  of  length-to-diameter  ratio  is  presented 
in  graphical  form  for  several  pipe  diameters  in  Figure  2,  and  in 
Figure  3  the  mixing  lengths  for  various  initial  Reynolds  numbers 
are  shown.  All  the  presented  data  are  for  a  fixed  mass  flow,  pressure 
and  entrance  centerline  temperature.  The  results  show  that  as  the 
Reynolds  number  increases  the  mixing  length  increases.  The  adjustments 
of  the  temperature  profile  from  a  radial  variation  to  a  uniform  distri 
bution  at  succeeding  axial  stations  are  also  presented  in  Figures  U 
through  9.  The  variations  of  the  exponential  decay  factor  "8^"  with 
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Reynolds  number  are  given  in  Table  I. 

The  engineer  is  often  required  to  design  mixing  tubes  vith 
varying  degrees  of  vail  heat  losses.  Thus,  these  calculations 
give  the  length  of  a  perfectly  insulated  mixing  tube  for  turbulent 
flov.  This  is  the  minimum  length  mixing  tube,  since  any  vail  losses 
vill  require  a  greater  length.  Calculations  for  non-perfect ly 
Insulated  pipes  vere  presented  in  Reference  8.  In  the  solution 
for  non-perf ectly  insulated  pipes,  non-uniform  entrance  temperature 
profiles  vere  not  considered. 
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APPENDIX  A 


The  general  solution  to  Equation  lb  is 

g  *  A  F(a,b,c,t)  +  Bt1”0  P(a-c+l  ,b-c+l  ,2-c  ,t )  ( A— 1 ) 


where 


a+b 


6 

7 


and  c  *  ■£ 
7 


For  the  boundary  condition  given  by  Equation  8  to  be  satisfied 
at  z  «  o,  B  must  be  zero.  Therefore,  Equation  15  now  reduces 
to 

g  *  A  F(a,b,c,z7)  *  A  F^a,^  "  a)  »  j*  (A-2) 

But  at  z  =  1 ,  g  must  be  finite;  hence,  the  function  F  at 
z  *  1  is 


F(a,b,c,l) 


r(c)T(c-a-b)  . 

r(c-a)r(c-b) 


( A— 3 ) 


and  if 

b 


—  and  z 
7 


1 
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the  function  F  becomes 


F(a,  |  -a,  f,  1) 
7  7 


r(f)r($-°-f+a) 
r(f -a)  r(a) 


GO 


r(va)r(a) 


( A-U  ) 


Thus,  for  the  particular  values  of  b  and  c  given  above,  the 
series  for  the  function  F(a,b,c,l)  diverges  unless  it  terminates. 
This  requires  that  "a"  must  be  zero  or  a  negative  integer,  vhich 
leads  to  the  following  admissible  values  of  "a": 

a  =  o,  -1,  -2,  - ,  -n  (A-5) 

yielding  the  corresponding  eigenvalues 
-  >*9  (-»)  (4  ♦  n) 


where  n  *  0,  1,  2,  3  .  .  .  . 


(A-6) 


The  resulting  polynomial  solutions  are  included  in  a  more 
general  class  known  as  Jacobi  polynomials  and  are  defined  by 
the  equation 


Fn(o;  y;  z7) 


F(-n,  1  +  n;  z7 )  m 

T  7 


( A-7 ) 
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I 


For  the  present  case 


a  *  n,  where  n  «  0,  1,  2,  3.  .  . 


a  *  - 


Y  *  — 


(A-8) 


Thus,  the  eigenfunctions  of  Equation  13  are  given  by  the 
equation 

A»  F-(fi  r  I?)  *  A» p”  (-“•  f +  ni  r  I?) 


(A-9) 


where  the  are  arbitrary  constants. 


The  functions  Fn(z)  are  orthogonal  with  respect  to  z5  as  a 
weight  factor: 


s 

s 


1  z5  F  (z)  F„(z)  dz  =  o  for  m  n 
m  n 


(A-10) 


1  z5  F*(z)  dz  «  C*  when  m  *  n 
o 


In  terms  of  the  variable  t  *  z7,  the  orthogonality 
relation  is 


fl  tY  1  (l-t)a  Y  FmFndt  =  0  for  m  B  n» 
J  O 


( A-ll ) 
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where 


and 

I 


above 

s 


where 


Y 


6 

7 


a 


6 

7 


1  tY“1(l-t)a‘Y  F2dt 


r ( y )  r  (o  +  i-y) 
r(o) 


(a  +  1  -y)n 


n ! 


2n 


( A-12 ) 


It  is  easily  verified  that  for  the  values  o  and  y  given 
,  Equation  26  becomes 


1 

o 


z5 


F2dz  * 


(n!)2 


( A-13 ) 


■  a  (a  ♦  1)  (a  ♦  2)...(o  *  n-l)  for  n  1,  2,  3... 
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with  (a)  defined  as 
o 

(o)o  -  1. 

For  n  *  o,  1,  2,  3  and  4  Equation  27  yields 

C  *  0.4082 
o 

c  *  0.2609 

C  *  0.2155 
2 

C  *  0.1904 

3 

C  *  0.1738  ( A-l4 ) 

4 

The  first  five  Jacobi  polynomials  are  easily  obtained  by 
expanding  the  series 


F  (a;  y;  t)  ■  F(-n,  a  +  n;  y;  t)  *  1 
n 


n 


1(1} 


(a  ♦  n)  (  O  ♦  n  +  1) 


(g+n  +  k-  l)  * 


(A-15) 


k*l 


Y ( Y  +  l) . . . (y  +  k-1 ) 


for  n  *  0,  1  ♦  ...4  and  y  f  0t  -1,  ...-n  +  1. 
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Fop  n  *  o 


F  (a,  y,  t)  *  1 
o 

n  «  1 

F  (a,  y,  t)  «  1  -  2.1667  z7 
1 

n  -  2 

F  (a,  y.  t)  *  1  -  6.6667  z7  +  6.9231  z14  (A-l6) 

2 

n  ■  3 

F  (a,  y.  t)  «  1  -  13.5000  z7  +  35.3077  z14  -  2U.1269  z21 

3 

n  *  U 

F  (a,  y»  t)  *  1  -  22.6691  z7  +  107.2U29  z14  -  171.5886  z21 

4 

+  87.3830  z28 

With  the  Cn  and  Fn  thus  obtained,  we  now  proceed  to  obtain 
an  approximate  solution  to  Equation  7  by  applying  the  Galerkin 
technique. 

For  this  method,  it  is  convenient  to  define  a  set  of 

functions  ♦  (z)  such  that 

n 

FnU> 

♦  n(z)  =  - -  ■  n  -  1,  2,  ...  (A-17) 

C  r\ 

n  An 


21 


with  the  property  that 


n  J  o 


♦  («) 

m 


(A-i8) 


where 


Cm  o  if  m  j*  n 
[■  1  if  b  »  n 


For  n  *  o  we  shall  define 


(A-19) 


( A-20 ) 


while  the  eigenfunctions  *n(z)  for  n  *  1,  2,  3  and  U  are 
computed  by  substitution  of  the  An,  CQ  and  Fn  directly  into 
Equation  31.  Listed  below  are  the  first  five  eigenfunctions: 


♦  *  2.UU95 
0 

♦  *  0.U018  -  0.8706  z7 
1 

♦  «  0.2773  -  1.8U88  z7  +  1.9199  z14 

2  ►  ( A-21 ) 

♦  «  0.2206  -  2.9776  z7  +  7.7877  z14  -  5.321621 

3 

♦  «  0.1866  -  4.2293  z7  +  20.0080  z14  -  32.0128  z21 

4 

+  16.3028  Z28 
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FIG.  I  SCHEMATIC  DIAGRAM  SHOWING  THE  COORDINATES  USED 


Fig  2  CENTERLINE  TEMPERATURE  VARIATION  WITH  NON-DIMENSIONAL 
AXIAL  LENGTH 
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FIG.  6  RADIAL  TEMPERATURE  PROFILES  AT  DIFFERENT  AXIAL  POSITIONS 
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FIG.  8  RADIAL  TEMPERATURE  PROFILES  AT  DIFFERENT  AXIAL  POSITIONS 
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RADIAL  TEMPERATURE  PROFILES  AT  DIFFERENT  AXIAL  POSITIONS 


TABLE  I 


VARIATION  OF  THE  EXPONENTIAL  DECAY 
FACTOR  WITH  REYNOLDS  NUMBER 


R  (inches) 

Re  x  10  5 

B2 

n 

0.375 

2.2fc939 

1.63577 

5.04350 

10.20157 

20.37003 

0.500 

1.68704 

1.31831 

4.06469 

8.22172 

16.41676 

0.625 

1.34964 

1.11516 

3.43831 

6.95473 

13.88689 

0.750 

1.12470 

0.97263 

2.99888 

6.06589 

12.11209 

0.875 

0.96403 

0.86644 

2.67146 

5.40361 

IO.78969 

1.000 

0.84352 

0.78387 

2.41688 

4.88867 

9.76147 

1.125 

0.74980 

0.71760 

2.21254 

4.47534 

8.93615 

1.250 

0.67482 

0.66308 

2.04443 

4.13531 

8.25720 
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TABLE  I.  (Continued) 


R  (inches)  Re  x  10  5  B2 


1.375  0.613  0.61733 

1.90339 

3.85002 

7.68755 


1.500  0.562  0.57833 

1.78315 

3.60680 

7.20189 
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